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ABSTRACT.  Interferometry  with  the  NAVSTAR  Global  Positioning  System  (GPS) 
is  the  most  efficient  method  of  establishing  three-dimensional  geodetic 
control  on  local  and  regional  scales.  This  is  already  true,  even  though  the 
gpnstellation  of  satellites  is  incomplete.  We  consider  some  theoretical  and 
practical  aspects  of  using  differenced  carrier-phase  observations  of  the 
^S  satellites  to  establish  three-dimensional  control  networks.  We  present 
simple  and  efficient  algorithms  for  processing  multi-station,  multi-satellite 
observations.  These  algorithms  have  been  used  to  establish  three-dimensional 
control  networks.  Under  widely  varying  conditions,  we  have  obtained 
accuracies  of  1  to  1.5  parts  per  million  (ppm)  in  all  three  coordinates,  or 
3  cm  for  the  coordinate  differences  between  network  stations  separated 
by  10  to  20  km.  When  interferometry  is  used  also  to  determine  the  satellite 
orbits,  we  anticipate  improvement  of  the  accuracy  of  regional  control 
networks  to  the  level  of  0.1  ppm.  y  . 


security  classification  of  this  PAOEriWiF"  D»l»  Bnfnd) 


THREE-DIMENSIONAL  GEODETIC  CONTROL  BY  INTERFEROMETRY  WITH  GPS; 

PROCESSING  OF  GPS  PHASE  OBSERVABLES 

Y.  Bock,  R.I.  Abbot,  C.C.  Counsefman  III,  R.W.  King 
Dept.  Of  Earth  Atmospheric  and  Planetary  Sciences 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 

S.A.  Gourevitch 
Steinbrecher  Corporation 
185  New  Boston  St. 

Woburn,  MA  01801 

ABSTRACT.  Interferometry  with  the  NAVSTAR  Global  Positioning 
System  (GPS)  is  the  most  efficient  method  of  establishing 
three-dimensional  geodetic  control  on  local  and  regional 
scales.  This  is  already  true,  even  though  the  constellation 
of  satellites  is  incomplete.  We  consider  some  theoretical 
and  practical  aspects  of  using  differenced  carrier-phase 
observations  of  the  GPS  satellites  to  establish  three-dimensional 
control  networks.  We  present  simple  and  efficient  algorithms  for 
processing  mu  1 1 i -stat i on,  multi-satellite  observations.  These 
algorithms  have  been  used  to  establish  three-dimensional  control 
networks.  Under  widely  varying  conditions,  we  have  obtained 
accuracies  of  1  to  1.5  parts  per  million  (ppm)  in  all  three 
coordinates,  or  1  to  3  cm  for  the  coordinate  differences  between 
network  stations  separated  by  10  to  20  km.  When  i ntef erometry 
is  used  also  to  determine  the  satellite  orbits,  we  anticipate 
improvement  of  the  accuracy  of  regional  control  networks  to  the 
level  of  0.1  ppm. 


First  International  Symposium  on  Precise  Positioning 
with  the  Global  Positioning  System 
Rockv i I  I e,  Mary  I  and 
April  15-19,  1985 


INTRODUCTION 


GPS  interferometry  is  a  method  by  which  three-dimensional  relative-position 
vectors  between  observing  stations  can  be  estimated  with  respect  to  a  world-wide, 
crust-fixed  coordinate  system.  The  processing  of  GPS  observations  and  the 
accuracies  obtained  under  widely  varying  conditions,  not  only  for  relative¬ 
positioning  but  also  for  single-point  positioning  and  for  GPS  satellite  orbit 
determination,  have  been  described  by  Counseiman  et  al.  (1S83),  Bock  et  al.  <1984>, 
King  et  al.  <1384),  Bock  et  al.  <1965),  Abbot  et  al.  (these  proceedings),  and 
Ladd  et  al.  (these  proceedings). 


Me  assume  in  this  paper  that  the  basic  observable  is  a  precise  measurement  of 
the  difference  beween  the  phase  of  the  reconstructed  carrier  wave  of  the  signal 
received  from  a  GPS  satellite,  and  the  phase  of  a  reference  signal  generated 
within  the  receiver  at  one  site,  at  a  series  of  epochs  determined  by  the  clock  in 
the  receiver.  Me  call  this  observable  a  one-way  phase.  (Such  an  observation 
can  be  made  for  the  LI  and/or  L2  frequency  band.)  In  practice,  for  a  particular 
observing  session,  a  number  of  receivers  simultaneously  observe  several 
satellites  over  a  time  span  ranging  from  several  minutes  to  several  hours.  The 
simultaneity  of  observations  is  crucial  since  many  sources  of  error  tend  to  affect 
different  observations  equally  at  a  given  time,  particularly  instabilities  in  the 
receivers'  and  the  satellites'  oscillators,  and  errors  in  models  of  the 
satellite  orbits  and  the  propagation  media.  Differencing  of  the  simultaneous 
observations  at  different  stations  cancels  common-mode  errors.  Satellite- 
oscillator  errors  are  virtually  eliminated,  and  the  residual  errors  of  both 
atmospheric  and  orbital  origin  are  approximately  proportional  to  the  distance 
between  the  stations,  up  to  several  tens  of  kilometers.  We  call  the  resulting 
observables  sing! e-d i ff erences.  Differencing  again,  between  satellites  that 
were  observed  simultaneously,  results  in  an  observable  that  is  also  free  of 
receiver  oscillator  instability.  The  ability  to  observe  several  satellites 
simultaneously,  and  therefore  to  form  double-differences,  is  what  makes  GPS 
so  powerful  for  relative-position  determination.  This  fact  should  not  be  forgotten 
in  discussions  of  so-called  "undifferenced"  processing  methods. 


In  this  paper,  we  consider  some  theoretical  and  practical  problems  in  the  use  of 
multi-station,  doubly-differenced  phase,  observations  to  establish  three- 
dimensional  geodetic  control.  We  present  our  solutions  to  these  problems  in  the 
form  of  simple  and  efficient  algorithms.  A  companion  paper  describes  the 
application  of  these  ideas  to  GPS  satellite  orbit  determination  (Abbot,  et  al., 
these  proceedings). 


DATA  PROCESSING  ALGORITHMS 


Consider  a  typical  GPS  observing  session  in  which  M  stations  observe  up  to  N 
satellites  simultaneously.  Satellites  rise  and  set  throughout  the  observation  span. 
In  addition,  temporary  obstructions  of  the  satellite  signals  at  one  or  more  of  the 
tracking  sites  may  cause  the  number  of  one-way  phase  observations  to  differ  from 
site  to  site. 


Let  us  assume  for  a  moment  that  there  are  no  losses  of  data  so  that  matching 
sets  of  observations  are  collected  at  the  sites.  From  the  one-way  phases,  at  each 
epoch  one  can  form  M<M-l)/2  single-differences,  but  only  M-1  are  independent. 

The  single-differences  will  sum  to  zero  around  any  loop  of  stations,  regardless  of 
the  nature  of  the  observation  errors;  i.e.,  the  closure  is  trivial.  If  each 
baseline  were  estimated  independently,  the  baseline  vectors  would  sum  to  zero 
around  each  closed  figure. 


Differencing  the  one-way  phases  between  stations  introduces  50%  correlations 
between  the  single-differences  at  an  epoch,  in  the  simple  case  of  equa I -var i ance 
uncorrelated  errors  at  the  individual  sites.  If  these  correlations  are  accounted 
for,  any  choice  of  M-1  single  differences  will  yield  identical  baseline  estimates. 


With  this  point  in  mind,  consider  the  typical  observing  session  where  the 
observation  sets  differ  from  site  to  site.  The  following  two-step  algorithm  will 
extract  all  the  information  from  the  raw  observables  and  at  the  same  time  will  take 
into  account  the  correlations  that  appear  as  a  result  of  differencing  the  one-way 
phases. 


The  basic  (linearized)  observation  equations  for  one-way  phases  can  be  expressed 
in  familar  matrix  form  by 

(p  ~  AX+V+V^  (1) 

where  the  observation  vector  containing  the  differences  beween  the  observed 

phases  and  their  a  priori  computed  values,  A  is  the  design  matrix  containing  the 
partial  derivatives  of  the  phases  with  respect  to  the  parameters  of  interest, 

X  is  the  parameter  correction  vector,  V  is  the  vector  of  the  non-common-mode 
observation  errors,  and  V  is  the  vector  of  the  common-mode  observation  errors. 

For  the  estimation  model,  we  assume  for  the  (normalized)  non-common-mode 
observation  errors  that 
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(2) 

E(VV^}  =  I 

(3) 

whore  E  denotes  expectation  and  I  is  the  identity  matrix.  The  common-mode 
observation  errors  cancel  in  the  subsequent  differencing  and  can  be  neglected. 

The  assumption,  therefore,  is  that  the  non-common-mode  errors  of  the  one-way  phases 
are  statistically  independent  in  space  and  time,  with  equal  (and  normalized) 
errors.  The  parameters  of  interest  for  three-dimensional  control  applications  are 
discussed  in  the  next  section. 


Let  us  denote  the  one-way  phase  for  station  m  and  satellite  n  by(p(m,n). 

At  each  epoch,  M-1  uncorrelated  and  normalized  single-differences  can  be  computed 
for  each  of  the  N  satellites  by 


Dl(m,n) 


( m+ 1 ,  n  > 


1/2 


<4) 


using  a  Gram-Schm i dt  orthogona I i zat i on  scheme.  (Note  that  the  index  m  in  Dl(m,n) 
denotes  only  the  numbering  of  the  single-differences,  whereas  on  the  right-hand 
side  of  (4),  m  denotes  the  station  number.) 


For  example,  for  4  stations  observing  simultaneously  the  signals  from  satellite  1 


01(1,1) 


C0(2,l)-(^(1,1)] 
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01(2,1) 


01(3,1) 


C0(3,l) - (0(l,l)+(i(2,l))] 
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C(^(4,l) - ((^(1,1)+0(2,1)+0(3,1))] 
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(6) 


(7) 


It  can  be  easily  verified  that  under  the  assumptions  (2-3),  the  observation 
errors  of  this  set  of  single-differences  are  uncorrelated  (and  normalized). 


Suppos^that  at  some  epoch  the  view  of  station  1  to  satellite  1  was  obstructed 
so  that  ^(1,1)  was  unavailable.  Then  only  M-2=2  uncorrelated  single-differences 
could  be  formed  for  satellite  1  from  the  one-way  phases  collected  at  the  three 
remaining  sites  by 


01(1,1) 


01(2,1) 


C(^(3,l)-<^(2,1)3 
C(^(4,l) - (C^(2,l)+(^(3,1)): 
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(8) 


(9) 


This  simple  scheme  ensures  that  all  one-way  phase  observations  are  used,  except  in 
the  case  when  only  one  site  is  observing  at  an  epoch. 


However,  the  one-way  phases  contain  information  on  the  absolute  positions  of  the 
observing  sites,  as  well.  Therefore,  if  we  include  at  each  epoch  the  common-mode 
observable  for  each  satellite 
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then  the  coordinates  of  each  observing  station  can  be  estimated  with  respect  to  a 
global  terrestrial  reference  frame.  Of  course,  the  relative  positions  of  the  sites 
will  be  more  accurately  estimated  than  the  absolute  positions  <see  Bock  et  al., 
1984).  Note  that  the  common-mode  observable  is  orthogonal  to  the  single¬ 
differences.  For  receivers  with  crystal  oscillators,  a  between-sate I  I i tes 
observable  (see  below)  is  more  valuable  for  point-positioning  since  any 
instabilites  in  the  receivers'  oscillators  are  eliminated.  The  one-way  phase 
common-mode  observable  is  valuable  for  more  accurate  frequency  standards  (e.g., 
ces i um-beam) . 


For  the  most  accurate  baseline  determination,  a  second  difference  is  taken, 
between  satellites.  The  above  orthonorma I i zat i on  scheme  can  be  applied  again. 
The  double-differences  are  denoted  by  02<m,n)  where  n  is  an  index  that  runs  from 
1,  N-1,  where  N  is  the  number  of  satellites  observed  simultaneously  by  the  M 
stations.  Me  form  for  each  set  m  (see  (4))  of  single-differences 


(Note  that  the  index  n  in  D2(m,n)  denotes  only  the  numbering  of  the  double¬ 
differences,  whereas  on  the  right-hand  side  of  (11),  n  denotes  the  satellite 
number.)  Thus,  for  M  stations  observing  N  satellites,  there  are  (M-1)(N-1) 
orthonorma I i zed  double-differences  that  can  be  formed  per  epoch.  This  scheme 
is  applied  to  all  visible  satellites,  so  that  risings  and  settings  of  satellites 
are  handled  naturally  (as  obstructions  are  handled  in  the  first  orthonormalization 
described  above).  Note  that  the  order  of  orthonormalization  could  be  reversed. 
That  is,  the  satellite  orthonormalization  could  precede  the  station  orthonormal¬ 
ization.  The  latter  order  may  be  preferred  for  point-positioning  with  satellite- 
differenced  observations. 


PARAMETER  ESTIMATION 


The  mathematical  model  relating  the  one-way  phases  to  the  geodetic  parameters  can 
be  found  elsewhere  (e.g..  Fell,  1980;  Remondi,  1984).  More  we  point  out  that  the 
partial  derivatives  used  in  a  weighted  least-squares  fit  of  the  orthonorma I i zed 
double-differences  are  manipulated  in  the  same  way  as  the  observables  in  (4)  and 
(11)  (i.e.,  they  are  also  orthonorma I i zed  as  the  normal  equations  are  being  filled). 

In  three-dimensional  network  applications,  for  M  stations  observing  a  total  of 
N  satellites,  we  adjust  as  a  minimum  3(M-1)  station  coordinates  (3M  if  we  are  also 
point-positioning  using  one-way  or  satellite-differenced  phases)  and  (M-1)(N-1) 
phase-bias  parameters.  The  phase-bias  parameters  represent  the  ambiguities  of  the 
one-way  phase  observables.  If  the  receiver  clocks  were  not  synchronized  and  we 
do  not  know  the  departure  from  synchronization  a  priori  (e.g.,  from  field 
measurements),  then  we  solve  for  M-1  clock  offset  parameters  as  well.  Rate  and 


higher  order  clock  terms  cannot  be  determined  well  from  the  double-differences. 
They  are  better  determined  from  the  single-differences. 
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Theoretically,  each  of  the  double-difference  phase-bias  parameters  has  an  integer 
value.  (The  value  would  be  zero  but  for  the  intrinsic  ambiguity  of  the  phase 
observable  —  that  is,  the  i nd i st i ngu i shab i I i ty  of  the  cycles  of  a  periodic 
signal).  In  general,  the  most  accurate  baseline  estimates  are  obtained  if  these 
biases  can  be  constrained  to  the  correct  integer  values.  If  the  statistical 
significance,  or  level  of  confidence,  of  this  integer  determination  is  sufficient, 
we  adjust  again  with  the  bias  parameters  fixed  at  the  best-fitting  integer  values. 
We  call  this  a  b i ases-f i xed  solution.  Otherwise,  we  consider  the  preliminary 
adjustment  with  its  real-valued  biases  as  final.  We  call  this  a  b i ases-f ree 
solution. 


With  a  typical,  three-  to  five-houi — long,  schedule  of  observations  from  single¬ 
band  GPS  instruments,  and  with  the  present  0.5-  to  1-ppm  level  of  uncertainty  in 
the  satellite  orbits,  we  find  that  we  can  fix  biases  when  the  baseline  length  is 
up  to  20  or  30  km.  Strange  (1384)  has  reported  success  in  fixing  biases  on 
baselines  of  up  to  50  km  in  length.  With  improved  satellite  ephemerides  and/or 
with  the  use  of  dual-band  receiving  instruments,  fixing  biases  is  possible  for 
longer  baselines  and  for  shorter  observing  schedules. 


PRACTICAL  DIFFICULTIES  AND  SOLUTIONS 


A  practical  difficulty  with  double-orthonormalization  in  the  mu  1 1 i -stat i on  and 
mu  1 1 i -sate  I  I i te  mode  stems  from  the  occurrence  of  c vc I e-s I i os  in  the  one-way 
phase  measurements.  A  cycle-slip  is  an  occasional,  sudden,  gain  or  loss  of  some 
whole  number  of  cycles,  due,  for  example,  to  temporary  occultation  of  a  satellite. 
If  the  slip  is  not  corrected  in  the  data  processing  by  the  addition  of  the  correct 
integer  number  of  cycles  for  the  appropriate  station,  satellite,  epoch  (and 
frequency  band  in  the  case  of  dual-band  observations),  then  obviously  the  geodetic 
position  estimates  will  be  corrupted. 


The  ability  to  recognize  a  cycle-slip  and  to  be  certain  of  "fixing"  it  depends 
on  the  magnitude  of  the  "noise"  in  the  observations,  relative  to  the  magnitude  of 
one  cycle.  For  a  baseline  of  the  order  of  10  km,  for  example,  the  r.m.s  noise  leve 
in  the  double-difference  observations  is  typically  about  0.1  cycle.  Therefore  a 
slip  of  just  one  cycle  is  conspicuous  in  a  series  of  double-difference  observations 
from  a  single-pair  of  satellites  and  a  single-pair  of  stations.  But  if  obset — 
vations  of  multiple  satellites  from  multiple  stations  are  linearly  combined  in  a 
double-orthonormalization,  a  single-cycle  slip  in  one  of  the  original  observations 
becomes  mapped  into  f rac t i ona I -eye I e-s I i ps  in  many  new,  orthonormal  observations. 
With  many  satellites  and  many  stations,  the  number  of  combinations  is  enormous. 


Our  approach  to  the  problem  of  cycle-slips  is  to  introduce  automatically,  in  a 
preliminary  adjustment,  a  new  phase-bias  parameter  whenever  a  gap  is  encountered 
in  the  phase  data.  Most  cycle-slips  (particularly  the  larger  ones)  are  associated 
with  losses  of  data  (gaps)  at  one  or  more  epochs.  The  post-fit  residuals  of  the 
(simple,  not  orthonorma I i zed )  double-differences  indicate  the  exact  epoch. 


.  •  ft  . 


satellite  and  frequency  band  where  a  cycle-slip  occurred.  (This  preliminary 
adjustment  also  serves  to  provide  good  starting  estimates  for  the  baseline  vector 
components.)  Now  the  cycle-slips  are  fixed  appropriately  in  the  one-way  phase  data 
and  a  second  adjustment  is  performed  using  the  double-orthonormalization  algorithm 
on  the  cycle-slip-free  data.  The  integer  phase-biases  are  then  fixed,  if  possible, 
as  discussed  earlier. 


Note  that  it  is  necessary  to  double-difference  the  phase  data  (1)  in  order  to 
correct  the  one-way  phases  for  cycle-slips  (since  the  noise  level  is  sufficiently 
low  to  detect  the  cycle-slips  only  in  double-differences)  and  (2)  in  order  to 
fix  the  phase-biases  at  their  best-fitting  values  (since  the  biases  are  integers 
only  after  doub I e-d i f f erenc i ng ) . 

For  an  application  of  the  approach  outlined  in  this  paper,  see  the  paper  by 
Bock  et  ai.  (1385),  in  which  the  analysis  of  a  35-station  three-dimensional 
geodetic  network  is  described.  The  results  of  that  analysis  indicate  that  three- 
dimensional  geodetic  control  can  be  routinely  obtained  at  the  1-2  ppm  level  under 
widely  varying  conditions.  When  interferometry  is  used  also  to  determine  the 
satellite  orbits,  we  anticipate  improvement  of  the  accuracy  of  regional  geodetic 
control  to  the  level  of  0.1  ppm. 
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